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Abstract 

This study shows that isoscaUng, usually studied in nuclear reactions, is a phenomenon common 
to all cases of fair sampling. Exact expressions for the yield ratio R21 and approximate expressions 
for the isoscaling parameters a and (3 are obtained and compared to experimental results. It is 
concluded that nuclear isoscaling is bound to contain a component due to sampling and, thus, a 
words of caution is issued to those interested in extracting information about the nuclear equation 
of state from isoscaling. 

PACS numbers: PACS 24.10.Lx,02.70.Ns,24.60.-k,64.70.Fx,25.70.Pq, 25.70.Mn 



1 



I. INTRODUCTION 



Recent studies of the isospin dependence of nuclear reactions at intermediate energies 
have been studied by comparing fragmenting colhsions of similar mass and energies but 
different isospin. Experimentally, the ratio of isotope yields in two different systems, 1 and 
2, R2i{n,z) = Y2{n, z)/Yi{n, z), has been seen to follow, approximately, an exponential 
function of the neutron jumber, n, and the proton number, z, of the isotopes; dependence 
known as isoscaling 



-R21 = y2{n,z)/Yi{n,z) = Cexp(an + (1) 



where a and (3 are fitting parameters and C is a normalization constant. This fit has 
been studied under the light of several theoretical models of nuclear reactions and, under 
different assumptions (see eg. ^|), the fitting parameters can be expected to be related to 
the symmetry term of the equation of state, Cgym^ through 
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a 



^ [{Z,/A,f - {Z^/A.f] , (2) 

where T is the assumed temperature of both reactions Thus the interest on studying the 
isoscaling phenomenon: R21 has the potential of elucidating the behavior of Csym at varying 
isospin. Temperature, etc.. It is relevant to notice that according to the above displayed 
relation the coefficient a should approach as T increases. 

In a series of recent works, however, it has been shown that the phenomenon of isoscaling 
can be found very early -before thermalization- in classical molecular dynamics simulations 
of nuclear reactions [5], as well as in non-thermal physical phenomena, such as in percolat- 
ing networks. In particular, percolation in two ^^colors" , {i.e. protons and neutrons [6]), or 
in extended "polychromatic" nets has demonstrated that the isoscaling behavior of the 
form of equation ([1]) emerges as a direct consequence of simple combinatorial problems with, 
eg. a = \ia{q2/qi), where = Ni/Ai. These results point to the fact that isoscaling, although 
connected to the equation of state, can also be produced by non-thermal processes.and then, 
probabilistic aspects of the problem can play a disturbing role complicating the determina- 
tion of Csym from isoscaling deermination. Thus the motive of the present work: to study 
the probabilistic aspects of isoscaling. 

In this study we investigate the phenomenon of isoscaling in the very simplest scenario 
-free of the geometrical constraints imposed by bond percolation- of sampling "protons" 



and "neutrons" directly from an urn. In section [TTl we will obtain the correct isoscaling law 
using probabilistic arguments, and will show that equation ([1]) is a limiting approximation 
to the exact expression. In section IIIII we compare our results to experimental values, after 
which the manuscript closes listing several conclusions. 



II. ISOSCALING AND SAMPLING 



Consider the problem of building clusters containing a number a of nucleons by sim- 
ply grabbing these a particles from an urn in which there are A particles composed of Z 
"protons" and "neutrons" , i.e. A = N + Z. We assume the sampling to be without 
replacement, and with no interactions amongst the particles nor with the urn itself; these 
premises are known as simple random sampling {SRS) in the statistics circles. 

To use this setup to study isoscaling we first focus on determining the yield of fragments, 
Y{n, z), that the previous scheme would produce after a large number of samplings. Stands 
to reason that such yield would be directly related to the probabihty of drawing n neutrons 
and z protons, i.e. a = n + z. Repeating then for a second urn with a different isotopic 
composition, one can easily obtain the corresponding R21 and, thus, the scaling law. 

Let us first determine the probability of ending with a cluster composed by n neutrons 
and z protons, i.e. a = n + z. The number of ways in which a cluster of a particles can be 
obtained from randomly sampling an urn with A particles is (^) = A\/{a\{A - a)]). Out of 
these possibilities only x will correspond to clusters with n neutrons and z protons. 
Thus, if we now assume that the probability of getting a fragment of size a out of an urn 
with A particles is P(a, A), the normalized yield of such a sampling will be: 



Y{n, N, z, Z) = P{a, A)P{n, N, z, Z, A) = P{a, A) ^'' 2. ■ (3) 



In particular for the case of sampling a particles composed of n neutrons and z protons 
{a = n + z) from an urn with Ai = Zi + Ni the term P(n, A''i, z, Zi) becomes 

P{n,N„z,Z„A,) - - ^^pT^ipT^ A\ • 

And taking the ratio of this probability to the probability of obtaining n neutrons and 
z protons from a sampling of a nucleons from an urn with A2 nucleons composed of A''2 



3 



neutrons and Z2 protons, the isoscaling ratio is given exactly by 

P(a, A2)P(n,iV2,^,Z2,A2) 



R 



21 



P(a,Ai)P(n,iVi,^,Zi,Ai) 



(4) 
(5) 



which can be readily used for calculations taking P(a, A2)/P(a, A\) as an overall normaliza- 
tion. 

This expression, although close to the usual exponential law ([1]), is not a straight line in 
the linear-log plot of P21 versus A^; the fact that experimental data also deviates from such 
a linear behavior is reassuring (see eg. jsl). The approximate exponential law ([T]) can be 
obtained from the exact result (jl]) using the binomial approximation to the hypergeometric 
distribution (see eg. [9]) which, in this case, depends on the assumption that x <^ Xj for 
X = a, n, Xj = Aj, Xj, Zi and i = 1,2. The sampling isoscaling ratio then becomes 

where we have introduced the probabilities of extracting a neutron, 
proton. Pi = Zi/Ai. 

Comparing to equation ([T]) we identify the overall normalization constant as C = 
P{a, A2) / P{a, Ai) and the isoscaling parameters as a{S) = ln(g2/q'i) = ^yi{N2Ai / NtAq) , 
P{S) = ln(p2/pi) = 111(^2^1/2^1^2), in perfect agreement with the percolation results p, l7|. 
Reviewing the procedure leading to equation ([6]), it is clear that exponential law is a direct 
result of the sampling. 

[For completeness, although not relevant to the nuclear case, we note that a functionally 
similar result can be obtained for the case of sampling with replacement. In this case, 
fragments with a = n + z will appear with probability P(n, N, z, Z, A) = Cp^q'"' = Cp^{l — 
pY~^, where C is given by the normalization C^^ = [p"^^^ — q^^^) /{p ~ q)j with p and q 
were defined before. Using this for urns 1 and 2 leads to 
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(Pi-gi)(-g2+'+P2+') 

where C(a) = [{2p2 - 1) +^1^^)] / [(2pi - 1) (-^2^^ + Ps^^)] • Since C[a) is inde- 

pendent of n and z, equation ([7]) depends on these variables in a functionally similar manner 
as the isoscaling law ([1]).] 



III. COMPARISON TO EXPERIMENTAL ISOSCALING 



The energy dependence of the isoscahng parameters has aheady been explored in coUi- 
sions [5I; here we compare the isoscahng parameter a obtained from experiments to those 
from samphngs, and study their variation as a function of beam energy. 

In view of the previous result, namely the fact that isoscahng can be expected from the 
mere act of fragmenting a system, the question to answer now is: what fraction of the nuclear 
isoscahng is due to sampling? 

In what follows we will show to which extent combinatorial effects (aka symentropy effects. 



see 



12j), are relevant by comparing the value of a obtained from experimental calculations 
with the one resulting from the simple combinaorial analysis. It is worth realizing at this 
point that in the case of the combinatorial analysis we are disregarding all correlations , in 
particular those associated with energy terms and then correspond to the very high Tem- 
perature limit. (see for example j3]).We then expect that the contribution of combinatorial 
terms to be apparent in collisions at high energies. Here we present a direct comparison 
of our findings to two sets of experimental data -not to try to reproduce them but- to 
assess the relative magnitude of the sampling contribution to isoscahng, and to attempt to 
draw a baseline from which future experimental studies will be able to extract the nuclear 
contribution to this effect. 



The data used was obtained by Yennello et al. at the Cyclotron Institute of the Texas 
A&M University 

The values of a{E) used in this comparison are shown in the inset of figured] and corre- 
spond to the ratios of the yields Y {^^ Ar + ^^Fe) /Y{^^Ca + ^^Ni), Y {^^Fe + ^^Fe) /Y{^^Ni + 
^^Ni), Y{'^°Ar + ^^Ni)/Y{^^Ca + ^^Ni) and Y{^^Fe + ^^Ni)/Y{^^Ni + ^^Ni) at the energies 
shown. The main panel of figure [1] shows the ratio of a{E) to the corresponding parameter 
obtained from the sampling, a{S) = ln(A^2^i/^i^2); cf. equation ([6]). 

It is easy to see that at all energies a{E) > a{S), and that for large energies a{E) a{S). 
A second point of interest is the fact that upon division by a{S), the previously scattered 
four curves of a{E) values collapse into a single group.. 
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FIG. 1: Ratio of experimental a{E) to the sampling a{S) as a function of energy. 

IV. CONCLUSIONS 

The conclusions of this work are quite simple. First we note that isoscaling can be 
expected to appear in any system undergoing a fair sampling; nuclear fragmentation -being 
a type of sampling- is bound to exhibit this phenomenon. 

Second, the inherent correlations of nuclear systems are expected to have an effect on the 
fairness of the sampling modifying the yield ratio R21 and the isoscaling parameters a and 
/3, as was demonstrated by a direct comparison to experimental results. 

Finally, a word of caution is needed if one attempts to extract information about the 
nuclear equation of state from isoscaling. Given that it is now known that a sampling- 
related isoscaling is ever present, obtaining quantities such as Cgym from equation [2] is not 
straightforward. In principle, the experimental results should contain the isoscaling produced 
by sampling, and -in some regime- both R2i{E) and a{E) should tend to R2i{S) and a{S). 

The indication that a{E) a{S) at high energies indicates -perhaps- that at those 
energies nuclear binding is less important and the ^^sampling" is closer to that of a non- 
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interacting urn; at lower energies, however, the reaction has longer interaction times and the 
phase space available for the sampling becomes a complex function of the energy distancing 
itself from the fair sampling case of equations (jlj) and (jH])- This is consistent with the 
implications of equation ([2]) which, through its inverse dependence on the temperature, 
indicates that in the limit of high energies the expected contribution from the equation 
of state to the isoscaling coefficients should vanish: a ^ 0. Likewise, studies of bond 
pe.co,a«on o„ pol.c>„on,atic ........ Q feene.ated .h.ough the nuCea. la«ce n,oda, 

at temperature T) have shown that in the limit of high temperatures only the probabilistic 
term survives. 
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